Low energy effective gravitational equations on a Gauss-Bonnet brane 
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We present effective gravitational equations at low energies in a Z2-symmetric braneworld with 
the Gauss-Bonnet term. Our derivation is based on the geometrical projection approach, and we 
solve iteratively the bulk geometry using the gradient expansion scheme. Although the original 
field equations are quite complicated due to the presence of the Gauss-Bonnet term, our final result 
clearly has the form of the Einstein equations plus correction terms, which is simple enough to 
handle. As an application, we consider homogeneous and isotropic cosmology on the brane. We 
also comment on the holographic interpretation of bulk gravity in the Gauss-Bonnet braneworld. 
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INTRODUCTION 



Recent progress in string theory provides a revolution- 
ary picture of the universe, in which our four-dimensional 
world is in fact a "braneworld" described by the motion 
of a thin wall embedded in a higher dimensional space- 
time. Various aspects on braneworld models have been 
addressed so far I*], and in a number of related publi- 
cations the simplest Randall-Sundrum model 2] is em- 
ployed for the analysis of gravity on the brane. In the 
Randall-Sundrum model, five-dimensional vacuum Ein- 
stein gravity with a negative cosmological constant is as- 
sumed. However, there will be higher derivative correc- 
tions in string gravity and hence it is natural to consider 
braneworld models with such corrections in the bulk ac- 
tion. While higher derivative corrections generally induce 
unwanted ill behavior, the special combination of curva- 
ture tensors called the Gauss-Bonnet term is known to 
avoid pathology: the Lagrangian is ghost-free, leading 
to well controlled field equations both at the classical 
and quantum levels 0- In this paper we consider on 
the braneworld model described by the Einstein-Hilbert 
action plus the Gauss-Bonnet term. 

In the context of the Gauss-Bonnet braneworld, cos- 
mology |3, ^ and linearized gravity 0, Q has been dis- 
cussed so far in the literature. In the present paper we 
consider nonlinear gravity in the braneworld. In this 
direction, Maeda and Torii 8] have derived the effec- 
tive gravitational equations on the brane using the ge- 
ometrical projection approach. Though in the Randall- 
Sundrum model the geometrical projection is a highly 
successful approach to clarify the nature of brane grav- 
ity 0, Maeda and Torii's effective equations in the 
Gauss-Bonnet braneworld are too complicated to be 
physically transparent. This is because the extrinsic 
curvature on the brane cannot be expressed solely by 
the energy-momentum tensor on it in a formal way. To 
cope with this problem, we shall resort to an approxima- 
tion method. The approximation employed here is that 
the relevant length scales on the brane are much larger 
than the bulk curvature scale, focusing on gravitational 
physics at low energies. 



The rest of the paper is organized as follows. In the 
next section we describe the model and present the basic 
equations following Ref. In Sec. Ill, we solve the 
bulk equation using the gradient expansion scheme and 
derive the effective equations on the Gauss-Bonnet brane. 
As the simplest example, we discuss homogeneous and 
isotropic cosmology on the brane in Sec. IV. Section V is 
devoted to summary and discussion. 



II. BASIC EQUATIONS 

The action we consider is given by 

+ [ d^xV^iQ + C„^-cJ), (1) 



where TZ is the Ricci scalar of the five-dimensional metric 
gMN and A is the cosmological constant in the bulk, h^i, 
is the induced metric on the brane, a is the brane tension, 
and Cm is the Lagrangian for matter localized on the 
brane. The Gauss-Bonnet term Cqb is defined as 



Cgb = 7^2 ^ 47^A/A^7^ 
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(2) 



and the parameter a has dimension of (length)^. The 
surface term Q is the Gibbons-Hawking-Myers term de- 
fined by [13 

Q^2K + ^a{j-2G;:K/), (3) 

where K,l^ is the extrinsic curvature of the brane and K 
is its trace, G,," is the four dimensional Einstein tensor 
with respect to the induced metric hf^^, , and J is the trace 
of 

1 



(4) 



In this paper we shall derive effective gravitational 
equations on a brane by solving the bulk geometry. For 
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this purpose, we use the (4+l)-decomposition of the met- 
ric and write the five-dimensional field equations in the 
form of the evolution equation along the extra dimension 
and the constraint equations. In deriving our basic equa- 
tions, we will closely follow the geometrical projection 
approach by Maeda and Torii (henceforth MT) (See 
Appendix for an alternative approach based on the work 
by Deruelle and Katz 11].) 

We write the bulk metric in the form 



gMKfdx'^^ dx^ — dy^ + q^y{y^x)dx^dx'^ , 



(5) 



where y is the fifth coordinate. We may assume that the 
position of the brane is y = without loss of generality, so 
that the induced metric on the brane is /i^tiy(x) — qfiv{y — 
0, x). We also assume a Z2-symmetry across the brane. 
The traceless part of the evolution equation is given by 

K,"K: - \5;k£kA , (6) 



where Kj^ is the traceless part of the extrinsic curvature 
K/(y,x) = (l/2)g"''5,(z^„, 



(7) 



is the "electric" part of the five-dimensional Weyl tensor, 
and is a hypersurface normal, n" — {dy)°'. is the 
traceless part of := M^^"" where 

M^auP ■= Rfiai^piq] — K^yKafj + Kf^pKya- (8) 

and Rf^avplq] is the four-dimensional Riemann tensor of 
the metric q^i,{y,x). 

The traceless part of the MT effective equations on the 
brane is arranged as 

^ (m/ + + a [(,) V + (,) V + (c) V] 



-— ^/A, 



(9) 



where M is the trace of . In the above we defined 



3 - aM 
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6(3 -I- aM) 3 + aM °' ^ 



(10) 



^va.^ fj, "I ^a.uj3^ 
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M^pE^I^M, 



(11) 



and 



with 



-^Mq^[aqfJ]u, (13) 

A^p^a Df,K^^- D^Kf,^, (14) 

iV^ := g^'^^oM/s = - D,,K, (15) 

and is the covariant derivative with respect to the 
four-dimensional metric q^v 

The trace part of the MT effective equations on the 
brane leads to the Hamiltonian constraint: 



1 



(16) 



where M Af/ - R[q] - + K^^K/. 
and (c)-ff are the trace of 

-8Af^„A4" - ^q^, {7M^ - 24Af„0Af"^ 
+3Af„/3,,Af"'5''^) , (17) 



+ -q^.M^pE'^l^ + 3ME^,, 



(18) 



and 



+2A^„/3^Af"'', + AN^o.pN^'^P 

+Zq^, (^A„A^" - lAr^^^iV"'37^ . (19) 

We require that an exact anti-de Sitter bulk with the 
curvature radius i solves the five-dimensional field equa- 
tions. Then the Hamiltonian constraint gives rise to the 
relation 



(20) 



where we introduced a useful dimensionless quantity 



(21) 



The Codacci equation, or the momentum constraint, 
is given by 

+2a (3 J/ - J<5/ - 2P^„''^V)] = 0, (22) 
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where 



+Rhf_i[ahp]^. (23) 



From the junction condition we obtain the following 
relation among the parameters: 



1 



boundary condition at the brane as 



equatic 



The (generalized) Israel junction equations give the 
or 



1 



^ 1 - ^ = ^n'cj. 



3 V2 



(24) 



where T^i, is the energy-momentum tensor on the brane. 

As a remark, we note that our basic equation @ holds 
on any y =constant hypersurfaces because the junction 
condition is not used in the derivation. 



III. EFFECTIVE EQUATIONS ON THE BRANE 

We will use the gradient expansion scheme [13. ITsj in 
order to solve the bulk evolution equations, and then will 
derive the four-dimensional effective equations describing 
the brane geometry. The small expansion parameter here 
is the ratio of the bulk curvature scale to the brane in- 
trinsic curvature scale, 



(32) 



B. First order equations 

The zeroth order result in the previous subsection is 
rather trivial; going to higher order we will obtain the 
field equations that govern the brane geometry hy^^{x). 

First let us look at the trace part of the extrinsic curva- 
ture, which can be determined without solving the bulk 
geometry. It is easy to see that 



{a) p 



\l>H^[l]H^O, (33) 



and then the Hamiltonian constraint reduces to a simple 
form^ : 



(1-/3)(i)m = 0. 



(34) 



Thus we arrive at the same conclusion as the Randall- 
Sundrum braneworld: 



e^i^\R\. (25) 
For any tensor we expand it as 

a;^(-''>a; + ^'^a; + ('^a; + -- - , (26) 

where ^"^^ A^ = 0(e"). 

A. Zeroth order equations 

The zeroth order solution is given by 

V = -7V- (27) and 

Since we have 

^'^M^^-P = l^(-5;'Sf + SfSJ^), (28) 

(0)^"^ = 0, (29) 

one can confirm that Eq. (|27|l trivially solves the evo- 
lution equation, the traceless part of the MT effective 
equations, and the Codacci equation at zeroth order. It 
is easy to check that Eq. I|27|l satisfies the Hamiltonian 
constraint as well. The bulk metric is written in the form 

^"''gMNdx^Ux'^ = dy^ + a^{y)h^^{x)dx'' dx"" , (30) 

where the warp factor is given by 



(35) 



Let us move on to the investigation of the traceless 
part. It is straightforward to show that 

i^)fj^ - 1(0) ^(1)77 CQ7^ 



(6)— M 2 



MiT'^ - n (38) 



(i)M/ - + 7^'^ V + 7^''^'^/. (39) 



Then the traceless part of the MT effective equations 
reads 



(40) 



Substituting this into the traceless part of the evolution 
equation, we obtain 



(41) 



a = e 



(31) ^ We assume /3 7^ 1 throughout the paper. 
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where R^^lq] is the traceless part of the Ricci tensor 
R^^[q] = R^[h]/a'^{y). Note that Eq. gT)) is again ex- 
actly the same as the first order evolution equation in 
the Randall-Sundrum model. The above equation can 
be integrated to give 

(1) = -2^VW + ^^'^x/(^), (42) 

where '•^^X^TC^) integration constant whose trace 

vanishes: — 0- To avoid confusion, we explicitly show 
here that the Ricci tensor in the equation is constructed 
from the induced metric hf^^{x). Combining Eq. ^ and 
the trace part H35|) . we obtain 



i 1 f ^ 1 



The extrinsic curvature (|43|l is related to the energy- 
momentum tensor on the brane via the junction condi- 
tion. The junction condition at first order is 

Wk- + ^( -d) J - - (1) + -(Dp/ - i(i)p<5/ 
'3\2 £ £ 



2 3^ 



Eq. agrees with the result of at long distances. 
Note, however, that the field equations (|48|) are nonlin- 
ear. This indicates that nonlinearity of gravity does not 
affect the gravitational coupling at low energies/long dis- 
tances. 

The integration constant (^•'xjl is constrained from the 
Codacci equation: 



which implies that 



0. 



= 0, (50) 



(51) 



(44) 



Substituting Eq. (|43|) into the above we obtain 

Dj'h; = 0. (52) 

Namely, is transverse and traceless with respect 

to the brane geometry. This means that it corresponds 
to the dark radiation, and in fact we have ^^''E,'^ = 

~{2/£)'-^^Xa on the brane. In what follows we set 



Substituting 



(^-'Xu" = for simplicity. 



(45) 



C. Second order equations 



and 



(i)p ''°' = -iB''--S''R 



(46) 



into the junction condition, we can rewrite Eq. 14411 as 



We now go on to the second order calculations in order 
to see corrections to four-dimensional general relativity 
(other than ^^x/). 

First let us focus on the trace part. The Hamiltonian 
constraint implies that 



+Pi (r: - IsJ'r) . (47) 



6 



(1-/?)(2)M = -c 



1(1)^^2 _ 2(1) Af/^'^M/ 
6 ^ 



,(l)ra/3 (l)r P" 
^ per 



(53) 



Then Eq. (g7J) together with Eq. (021 yields the gravita- 
tional equations on the brane at first order: 



R'' -Is'^R^ - T " + -- ^(i)y " 

2^ e{l + f3) £l + (3 



Xu- (48) 



We can see that for vanishing (^'x/T Einstein gravity is 
reproduced at this order. 

From the coefficient of T,,'^ we can read off the four- 
dimensional gravitational constant at low energies as 

^2 



87rG = 



£{1 + P)- 



(49) 



In the linear perturbation analysis of Ref. ^ the scale- 
dependent gravitational coupling was obtained, and 



From the definition of M and L f we have 



and 



(1) r per _ ^ per 
^aS3 — '-'q/3 ' 



(54) 



(55) 



where Cappa is the four-dimensional Weyl tensor. Solv- 
ing for (^^if , we arrive at 



(2)^ 

J3_£^ 
1 - /? 24 



6 



(1)^2 _ (1)^ u(l)j^ p 



w - 2(i)if/(iUl/ + -^^Hi 

p " Q 



, (56) 
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where W is the trace of 



■Mil' — C f^^C 



(57) 



RecaUing that (^Uf = and ^^^M/ = i^/iY^hil' which 
is assumed to vanish, we obtain 



f f 1 

24 V ' 3 y 1 - /? 24 



W. (58) 



Next let us consider the traceless part. Again we need 
to solve the evolution equation to obtain the correspond- 
ing part of the brane extrinsic curvature. 

A simple calculation leads to 



(2) 



(a) 



(2) 



H.. 



(2)fV ly 



2W/ - 2(i)M/C^/« 

(i)M "(^^ M " --S "^^^M '^^^'^Mo° 

^2 1 - /3 ' 

M ^ M M 



0, 



(59) 

(60) 
(61) 



where we defined the traceless part of . Then 
the traceless part of the MT effective equations reduces 
to 



(2)^ 



_(2)^. 



3(1-/3) a' 



(62) 



Note here that we have l/a^{y) factored out so that 

is computed from the induced metric hfj_^{x). We again 

omitted the dark radiation term '^^^Xn ^^^^u ■ 



From the definition of M,^ we have 



(2)m'^ - (2) 6 !^ I (1) Z^- «(1) f 



(63) 



Using Eqs. and JHSJi we can write the evolution 

equation as 



£„(2)^^. = (2)^^. + ^(2)^^. _ (1)^(1)^^. 



/3£2 



2(1-/3) a' 



4 ' 



(64) 



The Ricci tensor (2)^^'' in Eq. ^ can be calculated as 
follows. Integrating the first order result l|43(l . we obtain 
the four-dimensional part of the bulk metric as 



a^h^^ + — (a2 - 1) 



,(65) 



where the boundary condition g^j^^^O, x) = is under- 
stood. Then the Ricci tensor is expressed as 



£2 



, (66) 



where := DaD" and one should notice that 5^" sat- 
isfies S,,'^ — and 



D.s;' = 0. 

Using Eq. (jfififl and the first order result 

£2 1 



^^^,R[h]R;[% 



(67) 



(68) 



we can integrate the evolution equation (|64|) . yielding 



(2)/?- 



v_ 

2 Va4 - 

/3^2 



£ 1 \ 1 ^3 - 



2(1 - /3) a- 



(69) 



where is an integration constant dependent only on 

the brane coordinates x^. To make the physical meaning 
of this integration constant clear, we define 



R^ R^ 



6(1-/3) ^ 

—RR.'f 
3 ^ 



(70) 



Now we can see that ^'^'^E^ = — (2/^)(2)t^^'^ on the brane. 
Thus we finally have the extrinsic curvature of the brane 
in terms of x^'': 



'SR^ R^ 



/3£3 



6(1-/3) 
which can be rearranged into 



W 5 

16^ 



(2)^ - 



(71) 



(2)X-_A-(2)J^ 



£3 



2 



/3£3 



6(1-/3) 



p 16 ^ 



+ ^'^X;.^ (72) 
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where 



+ls;R£Rf " ^V^'- (73) 

The form of V^^" is suggestive. Rewriting this in terms 
of a new variable defined by s^^ := — {1/2)5^ R, we 
obtain the same expression as the weh-known quadratic 
energy-momentum term which was first introduced in 
Ref. 13. 

The Codacci equation at second order is summarized 

as 



where 



(1-/3) {i^)K;-5;^)K) + ^y;^ 



0, (74) 



y^L ■- R-iia Rp ~ RR^L + 2^/^ 



24 



:5,rR' 



C^^^PRp--2V;. 



(75) 



The Codacci equation is formally integrated to give 



(1-/3) + 
= r/(x)+ci5/ + C2Z/ 



(76) 



where is the part that cannot be expressed in terms 
of local quantities. As for the local part, is the di- 
vergence free tensor defined by 

- U^R^ - D^D^R + SJ^D^R, (77) 



and ci and C2 are constant coefficients. Note that 
and are the two linearly independent, divergence 
free combinations of curvature tensors which are of or- 
der e^. This is because the variation of ^J—hR^ and 
^/ —hR^j.vR^" with respect to the metric gives Z^ and 
+ Z^ respectively. Due to the Gauss-Bonnet the- 
orem in four dimensions, another curvature squared term 
V—hRf_iupcrR^'^'"^ does not give rise to a new independent 
local tensor quantity. 

Substituting Eq. (O into Eq. (|7S|l . we obtain 



The traceless condition ^'^'^x^ = yields the constraint 



r/ = (1-3/3)- 

/3£3 



W - 3c2D^R. 



(79) 



In the cases of /3 = (the Randall-Sundrum model), 
T^'^ represents the trace anomaly, which is expected from 
the AdS/CFT correspondence in the braneworld 
model ,17|. The nonlocal nature of r^j^ indicates that it 
corresponds to the energy-momentum of the holographic 
CFT on the brane. However, it is not trivial whether 
this is the case when /3 7^ 0, because the AdS/CFT cor- 
respondence in the presence of the Gauss-Bonnet term 
has not been addressed much yet |T^ . 

Including the second order computation in the junc- 
tion condition, we shall derive the effective equations on 
the branc with corrections to four-dimensional general 
relativity. The junction condition now becomes 



(80) 



Using the results obtained in the previous and present 
subsections, we have 



-T, 



+ ^ ^(i + /3) 



— - (r/ + c,S; + c,Z;) , 



(81) 



G„ 



rp V 



21^(2).. 



3 



T_ 

16^ 



P 



(82) 



Here the constants ci and C2 are determined by the 
boundary condition other than that imposed at the 
brane. We stress that the above equations are correct 
even for nonlinear gravity, as long as they are applied to 
the low energy regime or at long distances. 



Ply. 

2 ^ 



+ (1-/3) 

pe^ 



(2; 



f ^ — -p ^ 



6(1-/3) 



W " S 



(78) 



IV. APPLICATION: COSMOLOGY ON THE 
BRANE 

As an application, let us consider homogeneous and 
isotropic cosmology on the brane. We write the induced 
metric as 



h^ydx^dx" = -dt^ + a^{t)5.ijdx''dx' , 



(83) 



This equation relates the integration constant '^'^^Xu to where a(t) is the scale factor. For this metric the Weyl 



the nonlocal part and the free parameters ci and C2- 



tensor vanishes, Cfiupa 



0, and hence we have W^"^ = 



7 



and = —2V^. Non- vanishing components of are 
given by 



tr4 



(84) 
(85) 



where H :— a /a and a dot stands for a derivative with 
respect to t. Then a straightforward calculation shows 
that 



D^y;- = -2D,v; = 0, 



(86) 



and thus the Codacci equation at second order implies 



(87) 



This means that '^-'x^," behaves as a conserved radiation 
like component on the cosmological brane (as is the case 
at first order). Assuming that this dark radiation term 
vanishes, we obtain the modified Friedmann equation 
from the [tt) component of the effective equations ()82|l : 



(1 - 3/3)^2 



3e{l+l3)' 4(l + /3) 



An exact form of the modified Friedmann equation in 
the Gauss-Bonnet braneworld was derived in Refs. 0,01 
which, omitting the dark radiation, is summarized as 



K^{p + (t) = 2\/H^ +e-^{3-f3 + 2f3fH^) . (89) 

Expanding the right hand side of Eq. (|89f) for small e — 
£^H^, we obtain the same expression as Eq. (|88|l . Thus 
the validity of the gradient expansion is confirmed. 

As long as one considers the homogeneous and 
isotropic universe on the brane, there is less advantage 
of our effective equations compared to However, it is 
sure that our procedure has a great advantage for general 
cases without symmetry. This is because we have the low 
energy effective theory at the nonlinear level. 



V. SUMMARY AND DISCUSSION 

In this paper we have obtained the effective gravita- 
tional equations at low energies in the "Gauss-Bonnet" 
braneworld. The derivation here is along the geomet- 
rical projection approach of Maeda and Torii At 
low energies we can solve iteratively the evolution equa- 
tion in the bulk by expanding the relevant equations in 
terms of e {^/L)'^ <C 1, where £ is the bulk curvature 
scale and L is the brane intrinsic curvature scale. Up to 
second order in the gradient expansion, nonlinear grav- 
ity on the Gauss-Bonnet brane is described by Eq. l|HT)) 
or Eq. H82I) , which is the main result of the present pa- 
per. Although the Gauss-Bonnet term makes the original 
governing equations complicated and far from transpar- 
ent, our effective equations have the form of the Einstein 



equations with correction terms, and are simple enough 
to handle. We can repeat the same procedure to include 
higher order effects. 

In the low energy effective equations (|81|l we have the 
nonlocal tensor t^i, and two free parameters correspond- 
ing to the bulk degrees of freedom. In the Randall- 
Sundrum model (/3 — 0), we can determine that part 
to close the equations via the AdS/CFT correspondence 
once we specify a CFT on the brane 0|. However, such 
a holographic interpretation of the bulk geometry is ob- 
scure in the presence of the Gauss-Bonnet term, and this 
point is an open question for further study. Not relying 
on the holographic argument, we can instead determine 
the integration constant Xnv by imposing another bound- 
ary condition, for example, at a second brane introduced 
away from the first brane. 
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APPENDIX A: COMMENT ON AN 
ALTERNATIVE FORMULATION 



The main text is based on the Gauss-Bonnet brane 
equations derived by Maeda and Torii 8] using a (4 4-1) 
decomposition first developed by Shiromizu et al. 9] to 
derive the Einstein brane equations. In this appendix 
we sketch a similar formulation of Gauss-Bonnet brane 
equations developed in Ref. in the case of Einstein 
gravity, with an emphasis on the structure of the equa- 
tions. 

Consider a five dimensional spacetime in gaussian nor- 
mal coordinates — {x^,y) with line element 



ds'^ = di/ + qf,^{xP ,y)dx^'dx'' 



(Al) 



and expand the metric coefficients q^j,y{x'' ^y) near the 
surface y = as 

+Oiy') . (A2) 
At lowest order in y the Riemann tensor of the met- 



8 



ric (jAip is 



where Dn and R 



2 Di^k^p 



1 



1 n k 



(A3) 



■fiupa 



are the covariant derivative and 



Riemann tensor of the metric g^^^x''). 

The Einstein-Gauss-Bonnet field equations are 



Expanding the Gauss-Bonnet tensor is only slightly more 
involved and yields 



A 



aUn = 



(A4) 



where Q^^ TZ^ — {1/2)5 j^TZ is the five-dimensional 
Einstein tensor and where the Gauss-Bonnet tensor 
is defined as 



'U A . n \-r,ALMN'r, 

rtg .— I K-BLMN 



1 



-27^^^7^BL + 7^7^/] - -Sb^Cgb- (as) 

Expanding the Einstein tensor at lowest order in y is 
a straightforward calculation which yields 



Gyy = [gp.f}9ua - .g^ii^ffa/?) k^"' k"''^ - ^i?, 

Gyp, = \ {9p.p9va - gpu9al3) D"" k"^'^ 



-D^ {k-p - S'^pk) 



(A6) 



(A7) 



Gp.i' — 7: i9pp9iya 9p}^9a[3) I ^ ~^ ^k^^ "rkkp^, 



(A8) 



Hy 



16 



Rpc.pR^'"'^ - ^Ro.pR'"' + R') (A9) 



a/3 



1 



\Dc. (^j/ - - 4P"^^,fc^^) , (AlO) 

-\N^.o.upr^ + 2M^, - ]^gp,C, (All) 



where P^'''P'' is defined in Eq. ^ in the main text, and 
the fully developed expressions of the other tensors are 



r 



-h'^pkp^k'^'' + h k^^pk''^ + h^'' {kp„k% - k^) , (A12) 

- (k^'^k""' ~ k'^'k'"^) + i (g'^'^g"" - 9'"'g°''^) {k^^k^p - fc^) _ k^'Pk^pg"" + k°'Pk^pgP'' 
-k-^Pk'^pgP''^ + kP-Pk^pg"!^ - kk°''^gP'' - kk^"" g°'f^ + kkP^g"!^ + kk'"^gP°' - 4PP°'''>^, (A13) 
3 3 1 ^ S 

-h°' i^T ka -J- —h'^ k-P h -I h"^ JcP h'^ k 4- —kh"' k^^ k a —hk" hP k 

^fi, ^jfi, ^h, ^rvfSj ^H. pt^ a-fv p)^va ^ pl^ cr'^fJ-f ^ g'^'^ fi'^ ^6 

k pkya + kpi/ + 2 (^RpaPjk'^j^k ^ Ruapjk'^pk ^) '^Rapk p^k^jj -I- —Rapk ^kp^ 

~~^Rpavpk'^^k^^ -|- —kRpaupk"'^ — k°'^ {RpakuP + Rwakpp) -\- —k (^Rpak°'^ -f Rvak°'^ 

+ -^Rp,(ik%kP, - k^) + ii? (2fc"^n - kkp,) , (AM) 



'gk"ijk'^^k'^^k\ + — (k'^pkP^) + kk"f^k^^k\ — -k^{k'^pkP^) + ^g^^ 
+Ro.f3jsk"^k''^ - AR^fsk^U^^ + 2kk'^^R^p + ^Rk%kP^ - ^Rk\ 

I 



(A15) 



There are no quadratic terms such as RpaR"^ in 'Hpi, as the Gauss-Bonnet tensor is identically zero in four 
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dimensions. As is well-known neither Tiyy nor Tiy^ 
contain Z^^; as for Ti^j^ it is quasi- linear in 1^,^. 

The Gauss-Bonnet field equations at lowest order in y 
can be seen as a boundary condition problem: given a 4- 
dimensional metric (7^,^ and an extrinsic curvature 
satisfying the constraints [the (yy) and (y/i) components 
of the field equations] , then the {fiv) component gives l^i, 
which is the metric in the bulk at order 0{y^). When 
a 7^ 0, Ifii, is only known implicitly because of the quasi- 
linearity of the Gauss-Bonnet tensor Hab ■ 

Now the junction conditions relate the extrinsic curva- 
ture to the energy-momentum tensor on the brane 
ll^ . They can in fact be read off from the (yii) compo- 
nent of the field equations, and are written in terms of 
k^u and as 




where the constant a := (6/f) [l — (4a/3^'^)] is such 
that fe^i^ = —{2/t)g^^, together with Rf^upa = 0, solves 
the brane equations l|A4ll . £ being related to A by A = 
-(6/^^) (1 - 2a/£2). (Here and hereafter we set = 1.) 
The (y/i) component of the brane field equations is then 
equivalent to the conservation of the energy-momentum 



tensor: -D^T^"" = 0. 

In Einstein gravity (a — 0), inverting the junction 
equations to express kp^u in terms of Tf^i, is elementary: 
k^iy = -(2/£).gp,y - [T^j^y - {l/'i)Tgf_,y\. Hence the brane 
field equations ljA4|l can readily be written in terms of 
T^n rather than /c^^and we arrive at the effective equa- 
tions on the brane [3 

ri — K^T — K — TP + —TT 

+ - Ya^^9pu. (A17) 

with the identification k\ — l/l. In deriving the above 
equation we introduced the electric part of the Weyl ten- 
sor and rewrote l^^, as 

= -2E^,y - ^Ag^^ -f ^fcppK^. (A18) 

In Einstein-Gauss-Bonnet gravity, the junction equa- 
tions (jA16p cannot explicitly be inverted to give k^y in 
terms of T^,j. Hence, the gradient expansion scheme 
developed in the main text, which consists in writing 
/cpi/ = —{2/l)g^y + e^i/ and expanding (jA16p to second 
order in e^j, , is a way to solve iteratively the brane equa- 
tions in a closed form. 
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